Divergent nematic susceptibility in an iron arsenide superconductor 
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Within the Landau paradigm of continuous phase transitions, ordered states of matter are char- 
acterized by a broken symmetry. Although the broken symmetry is usually evident, determining 
the driving force behind the phase transition is often a more subtle matter due to couphng between 
otherwise distinct order parameters. In this paper we show how measurement of the divergent ne- 
matic susceptibihty of an iron pnictide superconductor unambiguously distinguishes an electronic 
nematic phase transition from a simple ferroelastic distortion. These measurements also reveal 
an electronic nematic quantum phase transition at the composition with optimal superconducting 
transition temperature. 

PACS numbers: 74.25.F-, 74.25.fc, 74.25.N-, 74.70.Xa, 75.47.-m, 75.60.Nt 



An electronic nematic phase transition refers to a 
phase transition in which the electronic system self- 
organizes with orientational order without forming spa- 
tial periodic order^lj]. For crystalline systems, such a ne- 
matic transition breaks a discrete rotational symmetry 
of the crystal lattice without altering the existing trans- 
lational symmetry. Canonical examples include half fill- 
ing Quantum Hall states [a, the field induced metamag- 
netic state in Sr3Ru207 3] and most recently the hid- 
den order state in URu2Si2 4], in which cases the crystal 
lattice's fourfold symmetry remains almost unperturbed 
yet the electronic ground states exhibit a strong two- 
fold anisotropy. [5i] . Recently both high Tc cuprates[6|-|9[ 
and iron pnictides |lCll4l2l | have been proposed as candi- 
date platforms that might harbour an electronic nematic 
phase, which opens up exciting new possibilities related 
to the interplay of nematic order with high temperature 
superconductivity. However, one of the key doubts ac- 
companied by the experimental discoveries is that the 
crystal lattice of these two systems does not retain a four- 
fold symmetry. In particular, in iron pnictides there is 
an orthorhombic structural distortion accompanying the 
rapid increase of resistivity anisotropy, which puts the 
legitimacy of the term "electronic nematic" into ques- 
tion. Here we report measurements of the resistivity 
anisotropy of Ba(Fei_a;CO:E)2As2 induced by a tunable 
uni-axial strain, which exhibits a divergent behaviour 
as the system approaches the phase transition from the 
high temperature side. Our result explicitly shows that 
the structural phase transition in Ba(Fei_a;Coa;)2As2 is 
purely driven by the instability in the electronic part of 
the free energy, and furthermore reveals an electronic ne- 
matic quantum phase transition at the composition with 
optimal superconducting transition temperature. 



crystal samples of Ba(Fei_:ECoa;)2As2 to probe the ne- 
matic response. As shown in Fig. [IjA), by gluing the 
sample on the side wall of a piezostack, strains can be 
applied by the deformation of the piezo, which is con- 
trolled by an applied voltage(Vp)[13|. The strain (i.e. 
the fractional change of length along the current direc- 
tion, ep = AL/L) was monitored via a strain gauge glued 
on the back side of the piezo stack. Both ep and the frac- 
tional change of resistivity (77 = Ap/po, where po is the 
resistivity of the free standing sample before gluing on 
the piezo stack) were measured at constant temperature 
while the applied voltage was swept, as shown in Fig. 
[IjB). The voltage dependence of rj and ep shows hys- 
teretic behaviour due to the ferroelectric nature of the 
piezo materials, yet the two quantities exhibit a linear 
relationship without any hysteresis(Fig. [HB)). The neg- 
ative slope of 7/(ep) indicates that the resistivity is higher 
along the shorter bonding direction, consistent with pre- 
vious results [ni H, El- 

The amount of strain transmitted to the sample {es) 
can be assessed by gluing another strain gauge on the top 
surface of the crystal, shown schematically in the lower 
panel of Fig. UK A). The comparison of 65 and ep for a 
Ba(Feo.955Coo. 045)2 As2 sample is summarized in Fig[T](C 
and D). For applied voltages \Vp\ < 150T^, the strain is 
fully transmitted to the sample for temperatures below 
approximately 100 K for typical thickness crystals (less 
than 100 /im). For lower voltages, |lp| < 50^, the strain 
is fully transmitted to even higher temperatures (Figure 
1(D)). The maximum strain that can be applied (|e| < 
5 X 10""') is substantially less than the lattice distortion 
developed below the phase transition(10^^ ^ 10^"^), and, 
as we show below, the system is always in the regime of 
linear response. 



We apply a tuneable in-plane uniaxial strain to single 
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FIG. 1: (Color online) (A) Schematic diagram of a piezore- 
sistance measurement (upper panel) and of a strain gauge 
measurement (lower panel). Details about the configuration 
are described in the supporting online materials. (B) (Upper) 
The relative change of resistivity (t\ = l^p/ po) of a BaFe2As2 
sample and the strain measured by strain gauge on piezo 
(ep = IS.L/L) as a function of voltage at T = 140A". The 
strain and resistance were measured along the [1 1 0]t direc- 
tion of the crystal. (Lower) Same set of data but r) is plotted 
against ep . The red line is a linear fit to the data. (C)(D) 
The difference of the strain between zero applied voltage and 
(C) Vp = l^W and (D) Vp = 50V". Dash lines indicate the 
temperatures below which the strain is fully transmitted to 
the sample. For low voltage this temperature window spans 
well above Ts for all compositions studied. 



vides a direct measure of the electronic nematic order 
parameter. Specifically, in the ordered state the resistiv- 
ity anisotropy ip — {pb — Pa)/{pb + Pa) can be used to 
define a nematic order parameter. For a strained crystal 
in the tetragonal state pb and pa refer to the resistivity 
in directions parallel and perpendicular to the applied 
compressive stress. It can be easily shown that rj = tp 
if the increase in pb equals the decrease in pa, and that 
the two quantities are directly proportional even if this 
is not the case. The same is also true for the derivatives 
of these quantities such that drj/de oc dtp /de (supplemen- 
tary online text). 

Representative data showing the electronic nematicity 
(ry) as a function of strain (ep) for a BaFe2As2 sam- 
ple are shown in Fig. [D^A) at various temperatures 
above the structural transition temperature T, . Data 
were fit by a straight line in a small range of strain 



o 



T = 140K 








T=150K^ 








T=160K VjW 






- 


T=180K ^J^A 








T = 210K ^^^&. 








■ T = 240K -^ "^"W^fc 




■ 


BaFe^As^ 

i 1 1 1 1 




S 


f* 



Bp(x10') 



100 



sr 
T3 




300 



300 



FIG. 2: (A) Representative data for BaFe2As2 showing the 
relative change of resistivity (77 = Ap/po) as a function of 
strain (ep = AL/L) at several temperatures above Ts- The 
nematic response was obtained by a linear fit of the data near 
zero applied vohage ( -5 x 10"^ < ep{V) - ep(0) < 1 x 10"*, 
indicated by the vertical dashed lines). (B) Temperature de- 
pendence of the nematic response d-q/dep. Vertical line indi- 
cates the structural transition temperature Ts= 138K. Red 
line shows fit to mean field model, described in the main text. 



near zero applied voltage. As shown in Fig. [2lD), 
the quantity drj/de, which essentially measures the ne- 
matic response induced by a constant strain, diverges 
upon approaching Tg from above. This divergent behav- 
ior is reminiscent of the resistivity anisotropy observed 
above Tg for samples held in a mechanical clamp jll|. 
However, as we explain below, there is an important 
distinction between measurements made under condi- 
tion of constant stress (mechanical clamp), and constant 
strain(measurement of drj/de in the current set up.) 

From the thermodynamic point of view the stress and 
strain are conjugate variables, and the stress (here de- 
noted as h) is the externally controllable force, whereas 
strain is the response of a mechanical system. Intuitively 



it might be more reasonable to regard stress as a sym- 
metry breaking field. However from the electron nematic 
stand point, stress only couples indirectly to the nematic 
order parameter through strain. This relationship can be 
best understood by the following Ginzburg-Landau free 
energy: 

F^^^'+^-i^^ + ^e^ + ^e^-XiPe-he (1) 

Here ip represents the electronic nematic order parame- 
ter, measured by the resistivity as discussed above, e is 
the elastic strain, and h is its conjugate stress, a, 6, c 
and d are the coefficients of the two order parameters in 
the usual power series expansion, and A is the coupling 
constant. If there is a phase transition driven by the elec- 
tronic degree of freedom, then the coefficient a becomes 
zero at some temperature, i.e. a = ao(r — T*), whereas 
the other coefficients are temperature independent. The 
bilinear coupling term Xeip renormalizcs the coefficient of 
e^ (c becomes c - A^/(ao(T — T*))) such that the crystal 
lattice distorts simultaneously with the onset of nematic 
order. On the other hand if the phase transition is due 
to a structural instability, then it is the coefficient c that 
becomes zero {c = cq{T — T*)) [l^. Therefore the driving 
force can be distinguished by determining the tempera- 
ture dependence of the bare a and c coefficients. 

With this in mind we can now ask what is the difference 
between measuring the response of electronic nematicity 
Tp under constant strain e rather constant stress h. This 
can be answered explicitly by calculating the quantities 
of dip /dh and dip /de under the constraint of minimizing 
the free energy (supplementary online text): 
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From these expressions, it is clear that the nematic re- 
sponse under a constant stress (eq. ^ will show a 1/T 
divergence no matter whether the driving force is a struc- 
tural or electronic phase transition. However the nematic 
response under a constant strain will only diverge when 
it is a true electronic nematic phase transition(eq. [3]). 
In this sense, the divergence in drj/de shown in Fig. [2] 
is direct evidence that BaFe2As2 suffers a true electronic 
nematic instability, and the structural transition merely 
passively follows the nematic order. Since strain is a field 
to the nematic order parameter, we refer to the quantity 
dip / de as the nematic susceptibility. 

From Eq. ^dip/de = A/a = A/(ao(T - T*)), it is natu- 
ral to fit the data of drj/de in Fig, [2jB) with a Curie- Weiss 
temperature dependence. However Eq. [3] is only valid in 
the limit of vanishing strain, at which one can disregard 
the higher order non-linear terms. In the realistic exper- 
iment situation, there is always some intrinsic built in 
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FIG. 3: Temperature dependence of the nematic suscepti- 
bility of Ba(Fei_2,Co2,)2As2 for various compositions(open 
symbols). Successive data sets are offset vertically by 50 
for clarity. Solid lines are fits based on a phenomenologi- 
cal Ginzburg-Landau theory, taking into account an intrinsic 
built-in strain(supplementary online text) 



strain even at zero applied voltage due to the different 
thermal contraction of the sample and the piezo stack. 
To take into account this built in-strain, we perform a 
numerical fit based on the following expression: 
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~de^ aoiT - T*) + 'ihril 



Xo 



(4) 



The effect of next order non-linear term is included in 
the 3677q in the denominator, where 770 is the resistivity 
anisotropy induced by the built in strain as a function 
of temperature, measured by the difference of resistivity 
of a sample before and after gluing on the piezo stack. 
In addition to ag and h introduced before, xo is a fitting 
parameter to model the intrinsic piezoresistivity effect of 
the materials that is unrelated to the electron nematic 
phase transition. 

The result of this fitting is plotted in Figl^l^B) as a 
solid red curve, which is in excellent agreement with 
measured data drj/de. Interestingly, the mean field crit- 
ical temperature T* obtained from the fitting is 116ii', 
22K lower than the actual phase transition temperature 
{Tg^ 138K ). This can also be understood from the 
Ginzburg-Landau free energy in Eq. [TJ By minimizing 
the free energy it can be derived that the nonzero nematic 
and structural order parameters onset simultaneously at 
a temperature Ts = T* + AV(aoc), higher than T* . This 
is due to the bilinear coupling between the electronic ne- 
matic system and the crystal lattice, which lifts the crit- 
ical temperature of the electronic instability to a higher 
temperature. Physically, the lattice provides a polariz- 
able medium, which enhances the nematic instability. 

The divergence of drj/de not only reveals the ten- 
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FIG. 4: Evolution of the nematic susceptibility (drj/de) of 
Ba(Fei_xCoi:)2As2 as a function of temperature and dop- 
ing. Structural, magnetic and superconducting transition 
temperatures {Ts ,Tm and Tc) are shown as squares, trian- 
gles and circles. The mean field electronic nematic critical 
temperature(r*) obtained from the fit to the data in Fig. [3] 
are shown as open red stars. The evolution of nematic sus- 
ceptibility and nematic critical temperatures clearly indicates 
an electronic nematic quantum phase transition occurs close 
to optimal doping. 



dency towards an electronic nematic phase transition, 
but also measures the strength of nematic fluctuations, 
according to the fluctuation-dissipation theorem. We 
have measured drj/de of Ba(Fei_2:Co2;)2As2 samples for 
doping concentration ranging from the undoped par- 
ent compound to overdoped compositions, as shown in 
FigEl The magnitude of drj/de is plotted as a color- 
map in the composition verses temperature phase dia- 
gram in Fig|4l For the underdoped part of the phase 
diagram, drj/de increases rapidly near the structural 
phase transition boundary. As the doping concentra- 
tion increases, the intensity of fluctuations increases, and 
reaches a maximum near optimal doping concentration, 
where structural and magnetic transitions are fully sup- 
pressed. The nematic fluctuations persist to the over- 
doped regime, eventually decreasing as the superconduct- 
ing Tc decreases. The associated softening of the sheer 
modulus has been extensively studied by resonant ultra- 
sound measurements [12, llZl. 

To quantitatively track the evolution of nematic fluc- 
tuations across the phase diagram, numerical fits to the 
data were performed for each composition based on EqH) 
The obtained T* is also plotted as a function of compo- 
sition in Fig|4l It can be clearly seen that the mean field 
nematic critical temperature T* closely tracks the actual 
structural transition temperature Ts in the underdoped 
regime, and is suppressed to zero at the optimal doping. 
T* becomes negative as the doping further increases be- 



yond optimal doping, indicating a "paranematic" state. 
Significantly, our experimental data and analysis reveal 
an electronic nematic quantum phase transition for a 
composition close to optimal doping. It remains to be 
seen whether fluctuations associated with this quantum 
phase transition play an important role in enhancing T^ 
in the superconducting phase. Nevertheless the existence 
of nematic fluctuations across such a wide temperature 
and doping range suggests that they are a fundamental 
ingredient to describe the normal state of the system [12]. 

The experiment we introduced here is a methodology 
to detect the electronic tendency towards rotational sym- 
metry breaking as if there is no coupling to the lattice, 
and is by no means restricted to the iron pnictides. One 
can also incorporate this methodology with other exper- 
imental probes, which allows us to more generally disen- 
tangle the cause and effect in systems for which different 
degrees of freedom strongly couple. For example, it is still 
an ongoing debate about the microscopic mechanism of 
the electron nematic phase transition in pnicitides, which 
could potentially arise from the spin |18l422| or orbital 23- 
25[ degree of freedom. By applying a constant strain 



and measuring temperature dependence of the orbital 



response by ARPES[28l| or optical conductivity [26l |27| 
and measuring the spin response by neutron scattering , 
the debate might be resolved. The authors thank C-.C. 
Chen, P. Coleman, R. M. Fernandes, S. A. Kivelson, A. 
Mackenzie and Q. Si for helpful discussions. This work is 
supported by the DOE, Office of Basic Energy Sciences, 
under contract no. DE-AC02-76SF00515. 
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SUPPLEMENTARY ONLINE TEXT 
Materials and methods 

Crystal growth and transport measurements 

Single crystals of Ba(Fei_a;CO;i)2As2 were grown from 
a self flux, as described previously(Sl), and the cobalt 
concentration determined by microprobe analysis. The 
crystals have a plate-like morphology, with the c-axis 
perpendicular to the plane of the plates. The in-plane 
orientation was determined by x-ray diffraction. Electri- 
cal contacts were made using silver epoxy to sputtered 
gold pads, with typical contact resistances of 1 ^ 2 fJ. 
Resistivity measurements were made using a standard 
4-point configuration. 



Piezo device 

Motivated by previous piezoresistance measurements 
of quantum Hall systems(S2), we developed a similar 
scheme to apply in-situ tunable strain for single crys- 
tal samples of iron based superconductors. FigjS] shows 
a photograph of a representative commercial PZT piezo 
stack that was used for the experiment(S2). By applying 
a positive voltage bias through the red and black voltage 
leads, the piezo stack expands along the long dimension 
and contracts along the transverse direction. Samples 
were glued on the side wall using a standard two part 
epoxy (S3), and a strain gauge was glued on the other 
side to measure the amount of strain that is applied(S4). 
The relative orientation of the crystal axes with respect 




FIG. 5; Photograph of two representative crystals mounted on 
a piezo stack used to apply in situ tunable strain. Two sam- 
ples are mounted, corresponding to e// [110]t (upper crystal) 
and e// [100]t (lower crystal). Red/black wires are the pos- 
itive/negative voltage leads of the piezo. The (0 1) surface 
of the crystal is exposed, enabling transport measurements. 



to the piezo stack determines the direction in which the 
strain is applied. For a control experiment, two samples 
of Ba(Fei_a;Coa;)2As2 with x — 0.065 (i.e. nearly op- 
timally doped) were measured, one with strain aligned 
along the tetragonal [110] direction(e//[110]T) and one 
aligned along the tetragonal [100] direction(e//[100]T). 
The resistivity was always measured along the strain di- 
rection. As shown in Fig. [51 for e//[110]T the temper- 
ature dependence of drj/de exhibits a similar divergent 
behavior to that of the optimal doped sample shown in 
the main text FigH] . However, for e//[100]T (i.e. at 45 
degrees to the orthorhombic direction observed for un- 
derdoped compositions) drj/de has an opposite sign and 
a much smaller value, and exhibits a much weaker tem- 
perature dependence. The significant difference of drj/de 
obtained between for the two orientations confirms the 
following analysis in which the divergent behavior is due 
to the coupling to nematic fluctuations. 



Resistivity anisotropy and nematic order parameter 

In general, the electronic nematic order parameter can 
be expressed as a function of the resistivity anisotropy 
-0 = {pb — pa) / {pb + Pa) ■ Taking a Taylor expansion for 
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FIG. 6: Representative data showing the tempera- 
ture dependence of the nematic response dri/dep of 
Ba(Feo.935Coo.o65)2As2 samples, for the two configurations 
shown in fig. SI. Samples were prepared so that the strain 
and current were aligned along the [f 1 0]t (Fe-Fe bonding 
direction) and [1 0]t (Fe-As bonding direction), shown as 
black and red symbols respectively. Vertical line indicates 
the superconducting critical temperature T^— 25K of these 
optimal doped samples. Inset shows the relative change of 
resistivity (rj = Ap/po) as a function of strain (ep = AL/L) 
at T = 30K. Lines show linear fits. 



small values of ip^ the leading term is linear, such that 
■0 is directly proportional to the electronic nematic order 
parameter when the order parameter is small. As dis- 
cussed in the main text, we measure the relative change 
of the resistivity rj = Ap/po of a strained samples in the 
tetragonal state. In the absence of strain Pa — Pb = Po ■ 
When an infinitesimal amount of strain 6e is applied, the 
resultant pa and pb (here b represents the compressive 
strain direction) can be expressed as a linear function of 
6e. 



^ 



Pa = /5o(l - aSe) 
Pb = Po(l + a'Se) 

^—5e + OiSe^) 



(5) 
(6) 

(7) 



If a equals a', it's straight forward to show that ip — rj, 
otherwise another proportional factor needs to be intro- 
duced in the relationship between these two quantities. 
If only the leading order term in the expansion of ^ in 
terms of 5e is considered, we then get: 



r n ; a + a 

2a' ' 

r n , a' + a 

oe > V ^ tp ~ T] 

2a 



(8) 
(9) 



In principle, the resistivity anisotropy will also be af- 
fected by acoustic phonon fluctuations associated with 



the structural transition. However for a fc = transi- 
tion the wavelength of such fluctuations is comparable to 
the sample size, and therefore scattering from such fluc- 
tuations makes negligible contribution to the resistivity 
anisotropy except for extremely close to T5(S5). 



Formalism of phenomenological Ginzburg-Landau 
theory 

In the presence of a crystal lattice, an electronic 
nematic phase transition breaks the crystalline point 
group rotational symmetry. For example in the case of 
Ba(Fei_2;CO:r)2As2 ,the C4 symmetry is lowered to C2 
symmetry. In this particular case we can write down the 
system's free energy in terms of an expansion of an Ising 
order parameter ip; 



2^ 



b ,4 



(10) 



To minimize the free energy one can take the derivative 
of F with respect to ip and set it to zero: 



dF 

dip 



aip -I- feV" = 



(11) 



This equation allows a non-zero solution for a < and 
6 > 0. Following the usual procedure, we express a by 
ao(r — T*) and set 6 as a positive constant to model the 
onset of a finite order parameter at temperature T* . In 
the presence of electron-lattice coupling, the non-zero or- 
der parameter in the electronic degree of freedom will al- 
ways induce a finite order parameter in the lattice degree 
of freedom. In the case of C4 to C2 symmetry breaking 
the lattice order parameter is measured by a orthorhom- 
bic strain e = (a — 6)/(a + 6). The effect of electron-lattice 
coupling can be modelled by the following free energy ex- 
pression: 

F = %P^ + y^ + '-e' + ^/ A^e (12) 

The term \ijje indicates that electronic order parame- 
ter and lattice order parameters are bilinear ly coupled, 
which is the lowest order coupling allowed by symmetry. 
We take the derivative of free energy with respect to "ip 
and with respect to e, and set them to be zero. 
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dF 
Ik 



aip + bip^ -Xe = 



ce + de^ - Xip = 



(13) 



(14) 



These two equations determine the electronic and lattice 
order parameters as a function of temperature under the 
constraint of minimizing the free energy. Substituting for 



e in equation [14] using the expression of e = (a'0- 
obtained from eq. [131 we get: 



- b%l]'^)/\ To see how system responds to stress, we take the deriva- 



similarly we get 
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(15) 



(16) 



Equation [15] has the same form as equation 1111 but the 
coefficients are renormahzed due to couphng to the lat- 
tice. In essence, the order parameter becomes non-zero at 

.2 

a temperature Ts higher than T*, where Tg == T* + ^. 
We can also look at equation [T51 which describes the lat- 
tice part of the free energy. Due to the coupling to the 
electronic degree of freedom, the coefficient c, which is 
the "bare" elastic modulus, is also renormahzed to a new 
effective modulus Ce// = c — A^/a = c — A^/(ao(T — T*)). 
Even though c might always be positive, meaning that 
there is no instability in the lattice part of the free en- 
ergy, the effective modulus still goes to zero at Ts, indi- 
cating the onset of non-zero spontaneous strain e. No- 
tice that the above derivation is completely reciprocal 
between the electronic and lattice order parameters - an 
instability in the lattice degree of freedom will also lead 
to non-zero electron nematicity below a structural phase 
transition. Then the question arises, how do we know 
that the electronic anisotropy in Ba(Fei_j;CO;r)2As2 is 
due to a real electronic nematic phase transition, rather 
than a parasitic effect of a pure elastic phase transition? 
A Gedankenexperiment one can do is to put the electrons 
in an infinitely stiff lattice, i.e. tune the parameter c to 
infinity. In this scenario if there is really an electronic ne- 
matic phase transition, then the lattice will always stay 
C4 symmetric while the electronic system will sponta- 
neously break C4 symmetry at temperature T*, which 
can be measured by the resistivity anisotropy. Another 
approach is to directly measure the bare coefficients a, 
and determine whether the bare nematic coefficient re- 
ally becomes zero at some temperature T* . This mea- 
surement can be done by perturbing the system with a 
symmetry breaking stress, and measuring the quantity 
dtp /de. To obtain this result, wc introduce stress ft., the 
conjugate field to strain e, into the free energy: 



„ « ,9 



\r 



Xxpe — he 



(17) 



The stress breaks the dual symmetry between ip and e 
in free energy, because "0 does not couple directly to h. 
Minimizing this free energy with respect to ip and e we 
get a similar set of "equations of motion" for tjj and e: 



dF 

dip 



dF 



ce 



a^ + bip'^ - Xe 
+ be^ -Xij-h 







(18) 
(19) 



five of equation [18] with respect to h: 

a^ + 3bip^— - X— = 
dh dh dh 

From this we get: 

d^ d-ip/dh X 

de de/dh a + Zbip'^ 

If we only consider the zero stress limit h ^ 
and the above equation reduces to : 

dtjj dtp/dh X 
de de/dh a 



(20) 



(21) 



0, then V) ~ 0, 



(22) 



Therefore if there is truly an electronic nematic phase 
transition, one should observe a Curie- Weiss diverging 
behavior of the quantity dip / de with a critical tempera- 
ture T*. As described in the previous section, our mea- 
surement of drj/de is directly proportional to dtp / de. The 
experiment is equivalent to measuring the temperature 
dependence of the nematic order parameter at constant 
strain, as can be appreciated from eq. 1181 Above Ts we 
can neglect terms proportional to ip'^ , leading to V' = Xe/a 
and hence dip/de — A/a. Although it is highly unphys- 
ical, if we also allow a direct coupling between mechan- 
ical stress and the nematic order parameter and include 
a term ghtp in the free energy, we obtain: 



d'ip eg + A 
de a + Xg 



(23) 



The separation of c and a into numerator and denomina- 
tor still allows us to distinguish the electronic instability 
from structural, since dip/de diverges only when a goes 
to zero. Physically, this distinction is possible because 
dip / de is related to the ratio of two thermodynamic quan- 
tities. We refer dip/de as a nematic susceptibility because 
in the limit that (7 = 0, the strain e is the only field to 
the nematic order parameter. 



Numerical fits based on phenomenological 
Ginzburg-Landau tlieory 

The measured drj/de was fitted based on the Ginzburg- 
Landau theory described in the previous section. Since 
the difference of thermal contraction between sample and 
piezo stack will always induce some finite built in strain 
even when the bias voltage is zero, the assumption of zero 
stress limit does not hold in the real experiment. Here, 
data are fit using the full expression for drj/de taken at 
finite rj obtained from Eq. [21] 



drj 
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ao{T^T*)+3brj^ 



(24) 
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FIG. 7: Temperature dependence of the relative change of 
resistivity rjo induced by the intrinsic built in strain associated 
with the differential thermal contraction of the piezo stacks 
and the crystals glued on their surface . Data presented here 
together with drj/de presented in the main text Fig. 3 were 
used for the numerical fit based on Ginzburg-Landau theory. 



sistently explained by an enhancement of electron-lattice 
coupling A. The enhancement of A not only increases 
A/ao, but also could potentially affect the "background" 
piezoresistivity effect measured by xo- However the nu- 
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FIG. 8: Doping dependence of the parameters obtained from 
the numerical fit based on the Ginzburg-Landau theory. 



where 770 is the induced relative change of resistivity of 
the sample at zero applied voltage comparing to before 
it was glued on the piezo. It is a function of temperature 
and was simultaneously measured when we performed 
voltage sweep at fix temperature in order to acquire 77(e). 
Representative data corresponding to the data of drj/de 
presented in Fig. [3] in the main text are shown in Fig. 
[71 A numerical fit was performed based on the following 
expression: 



drj/dt 



Pi 



T-P2 



P^vi 



(25) 



where drj/de and ryo are functions of temperature T, and 
are simultaneously fitted by four parameters. Pi, P2 and 
P3 correspond to A/ao, T* and 35/ao in the Ginzburg- 
Landau theory respectively. P4 = xo is a new parameter 
we introduced to take into account the intrinsic piezore- 
sistive response of metals and semiconductors that occurs 
due to the induced orthorhombicity even in the absence 
of nematic fluctuations. The obtained fit parameters are 
shown in Fig. |4] and Fig. |8] as a function of cobalt con- 
centration X. The contribution due to the built in strain 
is usually of the order of a few Kelvin, which will only be 
relevant if T is very close to T* . As described in the main 
text, T* changes sign near optimal doping, which is con- 
sistent with an electronic nematic quantum phase tran- 
sition. Inspection of FiglH also reveals an enhancement 
of A/ao and xo near optimal doping, which can be con- 



merical fit also becomes less consistent close to optimal 
doping (Fig. [S]), which might be due to the effect of 
quantum critical fluctuations. To distinguish the effect 
of electron-lattice coupling and quantum critical fluctu- 
ations is beyond the scope of this paper. 

Finally, the parameter 35/ao remains mostly constant 
across underdoped and optimal doped compositions but 
becomes negative for some overdoped compositions. This 
might be due to the fact that both drj/de and 770 become 
small for overdoped compounds, which makes the fitting 
less constrained. On the other hand for the overdoped 
compounds there is no phase transition at finite temper- 
ature, and hence the sign of parameter b is less crucial, 
and higher order coefficients can be introduced to bound 
the free energy. 
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